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Understanding of the puzzling phenomenon of high temperature superconductivity requires re- 
liable spectroscopic information about temperature dependence of the bulk electronic density of 
states. Here I present a comprehensive analysis of T— evolution of bulk electronic spectra in Bi-2212 
obtained by Intrinsic Tunneling Spectroscopy on small mesa structures. Unambiguous spectroscopic 
information is obtained by obviation of self-heating problem and by improving the spectroscopic res- 
olution. The obtained data indicate that the superconducting transition maintains the mean-field 
character down to moderate underdoping, and is associated with an abrupt opening of the super- 
conducting gap, which is well described by the conventional BCS T— dependence. The mean-field 
critical temperature reaches maximum at the optimal doping and decreases with underdoping. Such 
behavior is inconsistent with theories assuming intimate connection between superconducting and 
antiferromagnetic spin gaps, and support proposals associating high temperature superconductivity 
with the presence of competing ground states and a quantum critical point near optimal doping. 
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I. INTRODUCTION 

How does high temperature superconductivity (HTSC) 
evolve with decreasing temperature and what happens 
at the superconducting transition? Where is the real 
critical temperature Tc? Does HTSC becomes stronger 
or weaker upon approaching the undoped antiferromag- 
netic state? These highly debated questions are crucial 
for understanding the puzzling HTSC phenomenon. The 
answers to all those questions could be obtained from the 
analysis of temperature dependence of the superconduct- 
ing energy gap A(T) in the quasiparticle (QP) density of 
states. 

So far the majority of spectroscopic studies on HTSC 
were made by surface sensitive techniques [H HI, S 
H) B 0|- However, obtaining reliable spectroscopic in- 
formation from surface spectroscopy on HTSC is im- 
mensely difficult: atomic-scale c— axis coherence length, 
rapid chemical deterioration, presence of the surface 
states [1] and inherently different doping state of the 
surface, may preclude determination of bulk electronic 
properties by surface sensitive techniques. Furthermore, 
for Bi2Sr2CaCu208+5 (Bi-2212) the surface spectroscopy 
probes the BiO rather than the superconducting Cu02 
plane [7] . All this urges the necessity of bulk spectroscopy 
of HTSC. 

Intrinsic Tunneling Spectroscopy (ITS) provides a 
unique opportunity to probe bulk electronic properties of 
HTSC. This relatively new technique utilizes weak inter- 
layer (c— axis) coupling in quasi two dimensional HTSC 
compounds, in which mobile charge carriers are confined 
in Cu02 planes separated by some blocking layer (e.g. 
SrO-2BiO-SrO in case of Bi-2212). This leads to forma- 
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tion of natural, atomic scale intrinsic tunnel junctions, 
and to appearance of the intrinsic Josephson effect at 
T < [i, 1 [13, [m, [H, [H, [11, M, m- ITS is well 
suited for clarification of questions highlighted above: it 
is a direct spectroscopic technique, is not sensitive to 
phase coherence, has very high resolution [neV achiev- 
able), is mechanically stable and thus perfectly suited 
for T— dependent studies of HTSC (unlike surface probe 
techniques), and, most importantly, probes bulk elec- 
tronic properties of HTSC. 

This work represents a comprehensive analysis of tem- 
perature dependence of the bulk energy gap obtained 
by ITS on small Bi-2212 mesa structures. Unambigu- 
ous A(T) is obtained by careful obviation and cancela- 
tion of self-heating. Improved resolution by means of 
T— differential ITS allows tracing the gap in the phase- 
incoherent state at r > Tc- It is observed that at all 
studied doping levels, the superconducting gap opens 
abruptly in a mean-field manner and is well described 
by the conventional BCS temperature dependence. The 
mean-field critical temperature T™^ decreases with un- 
derdoping, thus confronting speculations about persis- 
tence of superconductivity up to very high temperatures 
above in underdoped HTSC. In most underdoped crys- 
tals a remaining, weakly T— dependent pseudogap (PC, 
Apg) is observed at T™-^' < T < T* ~ 120 - IbOK. 
The pseudogap seems to form a combined gap with the 
BCS-like superconducting gap below TJ"-^, and thus com- 
pete with superconductivity. No signature of the PC is 
observed at the optimal doping. The obtained results 
are inconsistent with theories assuming intimate connec- 
tion between superconducting and antiferromagnetic spin 
gaps, and support proposals associating HTSC with the 
presence of competing ground states and a quantum crit- 
ical point near optimal doping. 

In conventional low temperature superconductors 
(LTSC) superconducting transition occurs as a result of 
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the second-order phase transition. It is associated with 
abrupt appearance of the order parameter, represented 
by the superconducting energy gap A(T < T™^) cx 

1 — T jT™^ ^ and by hnear growth of the upper critical 

field: i?c2 oc cx X-TjT^^ . The correlation 77^2 oc 
is fundamental, because Hci is inversely proportional to 
the square of the coherence length ~ Cooper pair size, 
which is inversely proportional to the pair coupling en- 
ergy A. All this is perfectly described by the mean-field 
BCS-Eliashberg theory of superconductivity [l3| • 

But how does HTSC emerge with decreasing tempera- 
ture, and what happens at Tc remain unclear. For over- 
doped cuprates, thermodynamic characteristics reveal 
unambiguous evidence of the second-order phase tran- 
sition at Tc lisl . Similarly, analysis of Nernst effect p^ . 
equilibrium ]20l| an d fluctuation [U [l^] magnetization 
and resistivity [23[ reveal vanishing of Hc2{T Tc) in 
a wide doping range. However, in underdoped cuprates, 
characterized by the persistence of the normal state pseu- 
dogap at T > Tc, the superconducting transition at 
Tc becomes obscured and both thermodynamic [3 and 
transport [2^1 properties become abnormal. According 
to some reports neither A nor Hc2 [131 vanish at 
Tc, although different interpretations of similar data are 
possible ^'211, and vanishing gap at T ^ Tc was also 
reported '[2, M, 21, i28,] . 

A related controversy exists about doping dependence 
of the coupling strength in HTSC. Although Tc and Hc2 
decrease with underdoping [2^, the gap measured 
by surface sensitive techniques was reported to grow 
[l|, Q . This has been taken as evidence for a continuously 
increasing superconducting coupling strength oc A/Tc 
upon approaching the antiferromagnetic state, assuming 
an intimate connection between the two states faoi. If 
true, this would indicate that HTSC has a magnetic ori- 
gin. However, other experiments reveal the existence of 
two distinct energy scales, of which one indeed increases 
with approaching the antiferromagnetic state, while the 
other follows Tc at all doping levels 0, [H, [H, [13, HI] . 

Does HTSC becomes stronger or weaker with ap- 
proaching the antiferromagnetic state? Again, the an- 
swer could be obtained by understanding, what happens 
at Tc- If the coupling strength is increasing with un- 
derdoping, then so does T™^- Some researchers assume 
that T™-^ may exceed the room temperature already at 
moderate underdoping |5i, i24, 31]. To cope with the ap- 
parent decrease with underdoping of transport T^'*"'"^, at 
which phase coherence is achieved in transport measure- 
ments, one has to assume the existence of an extended 
region Tp'"''" <T < T™-^ in which the amplitude of the 
superconducting order parameter is large but the phase 
coherence is destroyed by thermal fluctuations [sl] . 

The extent of the phase incoherent state T™f — 
TPhase ^ c^T^f ig described by the Ginzburg-Levanyuk 
parameter Gi [S^. For clean LTSC the fluctuation re- 
gion is very small because of very small T™f /Tp ~ 10"*, 
where Tp is the Fermi temperature [s^. Even for HTSC, 



T^^/Tf < 0.1. Therefore, expansion of the phase- 
incoherent state well below T^^ requires Gi ~ 1, which 
can be achieved only by decreasing the dimensionality of 
the system [SS^]. Thus one has to assume that supercon- 
ductivity is either one-dimensional (e.g. due to the pres- 
ence of stripes [soj ) or zero-dimensional (as in granular 
superconducting films [sj])- In this case, there would be 
no second order phase transition, nor significant ampli- 
tude fluctuations of the order parameter upon establish- 
ing of the phase coherence at T^^"'^'^, the superconducting 
gap would persist at TP^'''"' <T < T^f , and could be 
directly measured by tunneling spectroscopy. Therefore, 
the knowledge of A(r) close and above J'p'"^^'^ is crucial 
for understanding HTSC. 

An important clue to understanding temperature evo- 
lution of electronic states in HTSC was provided by re- 
cent surface photoemission experiments [3], which, un- 
like earlier works, showed that A does not have a simple 
d— wave momentum dependence, but is described by two 
distinct energy scales The anti-nodal gap has weak 
temperature dependence at T < Tc and turns into the 
pseudogap at T > Tc- Furthermore, it increases with 
underdoping and tends to merge with the insulating gap 
in the undoped antiferromagnetic state. On the other 
hand the gap in the nodal "Fermi arc" region must be 
associated with superconductivity because it follows Tc 
and vanishes close to Tc at all doping levels. Although 
those observations are consistent with several previous re- 
ports [13, [11] , the reliability of surface spectroscopy 
of HTSC is now under question, because it fails to re- 
veal the electron-like Fermi surface, uncovered by recent 
quantum oscillation experiments j35j . 

The paper is organized as follows. In Sec. II I describe 
samples and emphasize the dramatic difference between 
surface and bulk properties of Bi-2212 crystals. Sec. Ill 
contains basic characterization of temperature and size- 
dependencies of ITS on small mesas. In sec. IV main 
experimental results on temperature dependence of the 
bulk gap in Bi-2212 are presented. It is shown that self- 
heating is effectively obviated by decreasing mesa size 
and can be simply canceled out from experimental data. 
The observed results are discussed and summarized in 
sec.V. In Appendix-A, self-heating and non-equilibrium 
effects in intrinsic tunnel junctions are analyzed. In 
Appendix-B artifacts of in-plane resistance and limita- 
tions on the junction size are discussed. 



II. EXPERIMENTAL 

Interlayer tunneling occurs in various layered HTSC 
compounds, such as Bi-2212 [11, [lEJp, [M [IS [11 , Bi- 
2201 Tl-2212[40, 41] YBaCuOgfTand some others 
[ist . as well as in intercalated compounds [11]. ITS was 
also expanded to non-HTSC layered compounds [H, \^ . 

Observation of the intrinsic Josephson effect [1] at T < 
Tc provides the most clear evidence for interlayer tun- 
neling in strongly anisotropic layered HTSC. At present 
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TABLE I: Parameters of studied samples: A'^— number of IJJs; _ onset of the resistive transition; _ appearance of 

the measurable critical current and establishing of the c— axis phase coherence; TJ"^- mean-field critical temperature, obtained 
from extrapolation of A(T); Asg(O) the superconducting gap at T — > 0; Apc the pseudogap at T^^; Uta- the c— axis thermal 
activation barrier in the normal state; Jc(0)- the critical current density at T 0. All samples were made from the same batch 
of Bi2Sr2Cai_iYa;Cu208+5 single crystals. 
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all major fingerprints of the intrinsic Josephson effect 
were observed, including Fiske d HO, im and Shapiro 
[13, nil steps in Current- Voltage characteristics (IVC); 
the Josephson plasma resonance 47] ; thermal activation 
and macroscopic quantum tunneling [l^ from the 
Josephson washboard potential; and the flux quantiza- 
tion [§, [13, [m, [H, [13] ■ The latter experiments explicitly 
confirmed the correspondence between the stacking pe- 
riodicity of intrinsic Josephson junctions (IJJs) and the 
crystallographic unit cell of Bi-2212. 

Several techniques for preparation of IJJs have been 
developed, such as patterning mesa structures on top of 
single crystals (48| , 3D-sculpturing by Focused Ion Beam 
(FIB) [ll|,[3 and double-side fabrication 

Here all measurements were performed on small mesa 
structures because they are best suited for ITS: they 
have the best thermal anchoring and are less prone to 
self-heating and other artifacts, as described in the Ap- 
pendix. All measurements were made in the three-probe 
configuration, which is more robust toward artifacts, as 
discussed in Appendix-B. 



A. Samples 

To avoid variations caused by different crystal stoi- 
chiometry, single crystals from the same batch of Y- 
substituted Bi2Sr2Cai_:cY3;Cu208+5 (Bi(Y)-2212) were 
used in this work, with the onset of superconductivity 
at optimal doping at T°"-^^^ ~ 96K. Obtained results are 
not specific for this batch, see e.g. Ref . [H, [13, [13] ■ 

Several mesas of sizes from '--^7x7to~3x iixm? 
were patterned simultaneously on each crystal by pho- 
tolithography and Ar-ion milling. All mesas usually con- 
tained the same amount of IJJs N . The virgin crystals 
were slightly overdoped with ~ %1K. Ar-ion milling, 
provides very uniform and controlled etching, but is ac- 
companied by substantial heating at high-vacuum. This 
results in partial out-diffusion of Oxygen, so that mesas 
become underdoped. To reduce self-heating, some mesas 



were trimmed to sub-micron size by FIB. Details of mesa 
fabrication are described in Ref . [l3j . 

Table-I summarizes crystals used in this study. To save 
space and avoid repetitions, I show only a limited amount 
of raw data and instead address the reader to previous 
works listed in Table-I. All of the studied mesas exhibited 
good periodicity of QP branches in the IVC, indicating 
good uniformity of IJJs in the mesas (see the insets in 
Fig. 2a) and references in Table-I). 

The solid line in Fig. 1 a) shows a typical resistive 
transition for a near optimally doped mesa on SMa sam- 
ple. Here i?o is the zero bias resistance measured with 
small IfxA ac current. One can see two transitions 
(three branches) in Ro{T) : the major part of the mesa 
goes into the superconducting state at ~ 92K, while the 
final transition takes place at ~ 40K. To clarify their ori- 
gins, the IVC of this mesa at 4.7 K is shown in the inset 
of Fig. la). It is seen that the IVC is strongly hysteretic 
and consist of periodic QP branches due to sequential 
switching of IJJs between superconducting and resistive 
states. The hysteresis is typical for tunnel junctions and 
is caused by low damping and large specific capacitance 
of the junctions [s^lssj. which allows junctions to remain 
in the resistive state even below the critical current. Dis- 
tinct branches in Ro{T) in the superconducting state are 
originating from different sections of the hysteretic IVC, 
as indicated by arrows in Fig. 1 a). 

Experiments on small mesas allow measurement of the 
QP resistance at different bias, i?*^^(/), in the super- 
conducting state [lll48l,!49|, which is otherwise shunted 
by the supercurrent |5J]. Circles in Fig. 1 a) show the 
zero-bias QP resistance Rq^ obtained by extrapolation 
to / ^ of the last QP branch with all IJJ's in the re- 
sistive state. The QP resistance at different bias can be 
also measured explicitly by first pulsing a current above 
the critical current Ic and then ramping it down to a de- 
sired value j49j . Bias yields an additional parameter for 
intrinsic tunneling studies, which may render crucial for 
correct interpretation of the data ^49. ] . 
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FIG. 1: (Color online), a) Zero bias resistance Rq vs T for a 
mesa on a nearly optimally doped crystal SMa. Inset shows 
I — V characteristics at T = 4.77^ which clarifies the origin 
of certain parts of the resistive transition, b) The same data 
shown as a thermal activation plot Ro/T (in a logarithmic 
scale) vs. 1/T. It is seen that in the whole normal state 
region T > Tc, Ro is described by the Arrhenius law (dashed 
line) with a constant TA barrier Uta — 25meV. 



B. Surface versus bulk properties 

As seen from the inset in Fig. la), the critical current 
of the first junction ~ 20fiA is much smaller than that 
for the rest of the junctions Ic ^ 270^A. This is the top 
IJJ in the mesa, between the two outmost CuO planes. 
It is seen that the surface CuO plane is superconducting, 
but has a lower critical temperature ~ 40K. 

Noticeably, the critical current and the QP resistance 
of the second junction, formed by the second and third 
CuO planes below the surface, is practically the same as 
for the rest of the junctions. This unambiguously shows 
that suppression of superconductivity is solely the surface 
phenomenon and occur only in the top CuO plane, and 
that "bulk" behavior starts already from the second CuO 
layer below the surface. 

The reduction of surface T^' in our mesas is predom- 
inantly caused by chemical deterioration in atmosphere 
during a short period between cleavage of the Bi-2212 
crystal and deposition of the top Au protection layer. 
Such deterioration was studied in detail in Ref. [sst , where 
it was shown that T[ could be increased to < 80K, 
if cleavage and deposition are made quickly without 
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FIG. 2: (Color online), a.) I — V curves at different T for a 
small mesa on a near optimally doped SMa crystal. The kink 
and transition to ohmic resistance at the sum-gap voltage is 
clearly seen. Inset shows multiple quasiparticle branches due 
to one-by-one switching of intrinsic Josephson junctions into 
the resistive state, b) dl/dV vs. voltage per junction for 
the same mesa. Arrows indicate the following characteristic 
features: the sum-gap peak, a minor double-gap dip at T < 
Tc, and a crossing point aX T > Tc- 



breaking vacuum (replicating conditions for surface spec- 
troscopy of HTSC 0, i, i, i, ). This is still substan- 
tially lower than the bulk Tc- 

The remaining suppression of of the surface layer 
is often attributed to the proximity effect with the elec- 
trode. However, I would like to note that although the 
top CuO layer has a lower (in some samples less than 
20K) the second layer is perfectly "bulk" , i.e. not de- 
teriorated with respect to deeper laying layers. Thus, 
there is no detectable proximity effect between the first 
and the second CuO layers. This is natural because the 
transparency of the interlayer barrier is low (prerequisite 
of the tunnel junction) and the c— axis coherence length 
is sub-atomic (< 0.3 A). Exactly for the same reasons 
there should be no considerable proximity effect between 
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the electrode and the top CuO layer, because those are 
also separated by the blocking BiO layer. 

More likely, the suppressed reflects the fundamental 
difference between the surface and the bulk e.g., because 
the surface is lower doped than the bulk. It should be also 
noted that even in UHV conditions the chemical deteri- 
oration of un-protected Bi-2212 surface is non-negligible. 
Note that the time of deposition of a monolayer is only 
~ 20s at the residual pressure p = lO^^Torr and even in 
state of the art surface spectroscopy systems is at best a 
matter of few hours. 

In any case the observed abrupt transition from the 
surface to the bulk properties within just one atomic 
layer from the surface clearly indicates that it is not at all 
granted that the surface spectroscopy can uncover sub- 
stantial information about bulk electronic properties of 
HTSC. All this urges the necessity of bulk spectroscopy 
of HTSC, as emphasized in the Introduction. 
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III. INTRINSIC TUNNELING 
CHARACTERISTICS OF SMALL MESAS 

Figs. 2 a) and b) show the c— axis I — V and dl /dV 
characteristics of a small mesa on the near optimally 
doped SMa sample at different T . A pronounced kink 
m. I — V (peak in dl/dV), followed by the ohmic and al- 
most independent tunnel resistance is seen p6j . Such 
the IVC is typical for superconducting tunnel junctions, 
as demonstrated in Fig. 3 a), and is associated with the 
sum-gap singularity at Vsg — 2A/e, providing the basis 
for ITS and opening a possibility to study bulk electronic 
spectra of HTSC. 

Fig. 3 a) shows numerically simulated IVCs at 
T =4.9K and 40 K for a superconducting tunnel junc- 
tion with a gapless density of states at the Fermi level. 
The gaplessness was achieved by introducing an appro- 
priate depairing factor F = 2meV into the BCS density 
of states. Parameters were chosen to fit experimental 
data for the S42 sample, shown in Fig. 3 c). Clearly, A 
decreased from 33 meV at 4.9 K to 29 meV at 40 K. Fig.3 
b) shows the same calculated IVC's in which both current 
and voltage scales are normalized by the gap. It is seen 
that the curves merge, because at low T < Tc/2 both V— 
and /— scales oc A. A similar scaling is also observed for 
experimental data, shown in Fig. 3 d), confirming that 
the kink is related to the sum-gap singularity. Additional 
arguments can be found in the Appendix A-1. 

The sum-gap singularity is not the only gap-related 
feature in ITS. Further check for self-consistency of our 
interpretation can be obtained from analysis of addi- 
tional, more subtle gap-related features in dl/dV. The 
arrow in Fig. 2 b) indicates a small dip at 214^. It was 
studied in Ref. [38| and was attributed to enhancement 
of non-equilibrium effects at eV > 4A, when relaxation 
radiation of tunneled quasiparticles becomes sufficient for 
breaking Cooper pairs. 



FIG. 3: (Color online). Comparison of a) simulated and c) 
experimental IVCs for S42 at 4.9 K and 40 K; b) and d) 
demonstrate collapse of the I — V curves, scaled by A(T') and 
indicate that both V— and /— scales are determined by A(T). 



A. Self-heating free characteristics 

It is fair to say that ITS has became a spectro- 
scopic tool as a result of reduction of self-heating. No 
T— independent ohmic tunnel resistance, as in Fig. 2 
a), could be seen in earlier works [1, IH, [5^, deahng 
with large structures. Larger self-heating in such struc- 
tures leads to development of an acute thermal instabil- 
ity [H^l at voltages much smaller than VsgN. Another 
problem in large structures is associated with the limited 
amount of supercurrent that the CuO plane can carry 
(see Appendix-B). It is interesting to note that exactly 
the same problems were encountered at the early stage 
of experimental studies on LTSC tunnel junctions [ssj . 

In order to obtain unambiguous spectroscopic infor- 
mation, the problem of self-heating has to be carefully 
addressed. The temperature rise due to self-heating is 
given by a simple expression [ssj 



ST^PRThiT), 



(1) 



where P — IV is the dissipated power and Rrh is 
the effective thermal resistance of the mesa, which is 
T— dependent and, therefore, bias dependent ^5^. De- 
tailed analysis of self- heating in ITS, including numerical 
simulations of distortion of IVCs by self-heating can be 
found in Appendix- A. 

In recent years different ways of obviating self-heating 
in ITS were employed, such as pulse measurements (36| . 
miniaturization [50j, |59| and heat compensation 37 1. In 
Ref. 50] it was emphasized that miniaturization decreases 
self-heating at a given voltage per junction and pro- 
vides an unambiguous way for discrimination of artifacts 
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FIG. 4: (Color online), a) dl /dV{V) (in a semi-log scale) for 
two mesas with different area on a moderately underdoped 
S92 crystal. It is seen how heating and in-plane resistivity 
bend upwards curves at high bias and reduce the sum-gap 
peak voltage for the large mesa, b) dl /dV{V) curves at differ- 
ent T for the small mesa. In is seen that the curves maintain 
the linear V-shape (in the semi-log scale) in the whole sub- 
gap region when self-heating becomes negligible. A sudden 
crossover at T = Tc from tunneling-like with T— independent 
slope, to thermal activation behavior with T— dependent slope 
is clearly seen. Data from Ref . [49i] . 



of self-heating (size-dependent) from electronic spectra 
(material property, size-independent). 

In Figs. 4a) and 5a) dI/dV{V) characteristics of large 
and small mesas on samples S92 (moderately under- 
doped) and SMa (near optimally doped), respectively, 
are shown. It is clearly seen how self-heating distorts the 
last QP branch in large mesas. At low bias, when self- 
heating is negligible, the characteristics of both mesas 
are undistorted and the log dl/dV curves remain linear 
and parallel, i.e., simply scale with the mesa area. How- 
ever, at larger bias log dl/dV in the larger mesas starts to 
growth faster (super-linear) because the sub-gap conduc- 
tance increases with T (see Fig. 2b). The larger mesas 
also reach the peak earlier, i.e., at lower voltage, as a 
result of suppression of the gap by T. 

Figs. 4 and 5 b) show T— evolution of logdI/dV{V) 
characteristics for the smallest mesas. It is seen that the 
log dl/dV characteristics for those mesas remain linear 
in the whole sub-gap region V/N < Vgg, indicating that 




FIG. 5: (Color online), a) dI/dV{V) (in a semi-log scale) 
at different T for the same mesas on near optimally doped 
crystal SMa before and after FIB trimming. The peak voltage 
is reduced in the larger mesa due to self-heating, b) dl/dV 
curves at different T for the small mesa. The characteristic 
V-shape in the semi-log scale is observed at all T. As in Fig. 
4, the crossover from T— independent to T— dependent slope 
occurs at Tc. 



they are not distorted by self-heating up to Vsg . For both 
mesas this was explicitly proven in Ref.[5^: for SMa by 
size-independence of the peak voltage and for S92 by in- 
situ measurement of self-heating. 

Note that heating-free ITS characteristics for both 
crystals have a remarkably trivial V-shape in the semi- log 
scale, and that the slope of the curves for both samples 
experience an abrupt crossover from thermal-activation 
(TA) like 1/T- dependence &i T > to tunneling-like 
T— independent slope at T < Tc. As discussed in Ref. [49], 
this indicates opening of an additional quantum trans- 
port channel for Cooper pairs at T < Tc. 



B. Thermal-activation behavior in the normal state 

In Ref. di] it was shown that at T > Tc ITS character- 
istics exhibit TA behavior. Up to moderately high bias 
they are described by a simple expression: 
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probably not the gap in electronic spectrum of CuO lay- 
ers, but the property of the blocking BiO layer. Pos- 
sible "non-gap" orig ins of the c— axis TA barrier were 
discussed in Ref. 43|. Those include resonant tunneling 
through the impurity state in the blocking layer, inelas- 
tic tunneling with excitation of a molecular mode in the 
barrier, and Coulomb blocking of tunneling in the poorly 
conducting two-dimensional electron system. As already 
noted in Ref.[2^, the later bares a striking similarity 
with experimental V-shape characteristics. The Coulomb 
blocking depends entirely on the conductivity of the two- 
dimensional electron system, which would naturally ex- 
plain the increase of TA barrier with underdoping. 



Improving resolution by T 
spectroscopy 



differential 



FIG. 6: (Color online). Simulated TA characteristics from 
Eq.(|3]) with Uta = 2AmeV . They reproduce all major fea- 
tures of experimental curves at T > Tc in Figs. 4 and 5 
b) including the V-shape with the slope proportional to the 
reciprocal temperature and the crossing point at eV ~ 2Uta ■ 



dl , , 1 
^(T,F)oc-exp 



Uta 



cosh 



eV 



2kBT 



(2) 



with a constant TA barrier Uta- Indeed, the cosh term 
reproduces the rounded V-shape of \n[dl / dV]{V) curves 
with the slope that increases as ~ 1/T. The TA behavior 
at zero bias at T > Tc is demonstrated in Fig. lb). In 
order to expand the TA-analysis to higher bias, we have 
to carefully integrate TA-current through the junction: 

/+00 
tTA{E,V)piE)p{E + eV) (3) 
-OC 

[/(£;)(! - f{E + eV)) - f{E + eV){l - f{E))] dE. 

Here txA = ex-p{~U{E,V)/kBT) is the transition 
probability of the TA process, and U is the effective TA 
barrier U{E,V) = mm[{UTA - E - eV/2);0], p{E) are 
density of states in electrodes and /(£')— Fermi distribu- 
tions. 

Fig. 6 shows simulated dI/dV{V) characteristics (in 
the semi- log scale) obtained from Eq.Q with Uta — 24 
meV and for completely normal electrodes p{E) =const. 
It reproduces all characteristic features of experimen- 
tal data at T > Tc, including the crossing point at 
eV ~ 2/7ta and the inverted parabolic shape at high 
temperatures ksT > Uta [131 ■ 

Apparently, the c— axis TA barrier should be identi- 
fied with the phenomenon referred to as the large c— axis 
pseudogap in the previous literature. However, the amaz- 
ing success of the trivial TA model, with only one con- 
stant parameter Uta in the whole normal region T > Tc 
and without any momentum dependent gap in the den- 
sity of states, suggests that the c— axis pseudogap is most 



As seen from Figs. 4 and 5 b), the sum-gap peak 
in dl / dV is rapidly smearing out with approaching Tc. 
Remaining weak spectroscopic features can be traced 
in a standard way by studying higher derivatives, e.g. 
d'^I/dV'^ d^. However, they are obscured by the para- 
sitic TA-background. In Ref. [131 it was shown that the 
ITS resolution at T < Tc can be improved by subtract- 
ing the TA curve at T > Tc. But since TA is strongly 
T— dependent, such subtraction does not completely re- 
move the changing TA background. 

Here I suggest the following optimization for TA back- 
ground cancelation: First, consider a normalized dif- 
ference between two characteristics as in Figs. 4,5 
b) at nearby temperatures T2 > Ti: T(Ti,T2) = 
[ln{dI/dV{Ti))-ln{dI/dV{T2))]TiT2/{T2-Ti). Accord- 
ing to Eg.©, for pure TA F(Ti,T2) ~ Uta - eV/2, i.e., 
is approximately T— independent, thus allowing optimal 
cancelation of the TA background. Another important 
advantage of this T— differential scheme is that it em- 
phasizes any T— dependent spectroscopic feature. 

In Figs. 7 a) and 8 a) T— differential characteristics 
T(Ti,T2) are shown for moderately underdoped crystals 
S92 and S82, respectively. It is seen that substantially 
above Tc the curves collapse into a single universal curve, 
as expected for pure TA. At lower T, the sum-gap peak 
and the double-gap dip are clearly resolved up to Tc. In- 
terestingly, small deviations from the universal TA curve 
are seen even above Tc up to ^ 12>QK . To see them more 
clearly, in Figs. 7 b) and 8 b) the universal TA curves at 
high Ti 2 were subtracted from the T— differential char- 
acteristics. In such the plot the TA background is com- 
pletely removed and we can very clearly see the remaining 
spectroscopic features close and even above T^. 



IV. TEMPERATURE DEPENDENCE OF THE 
BULK ENERGY GAP 

Figs. 9-12 summarize temperature dependencies of 
ITS features for different samples. Data for the same 
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FIG. 8; (Color online). The same as in Fig. 7 for another 
underdoped S82 crystal. The collapse of curves at high T into 
the universal TA curve is clearly seen in a). From panel b) 
it is seen that the normal state pseudogap remains roughly 
T-independent at T > Tc. 



FIG. 7: (Color online). Improvement of spectroscopic reso- 
lution at elevated T by means of T— differential spectroscopy 
for the same small mesa on S92 crystal. In a) shown are 
differences between ln{a = dl/dV) curves from Fig. 4 b) at 
two consecutive T, normalized by the temperature difference. 
Such curves emphasize T— dependent spectroscopic features 
and collapse into the universal curve for the case of TA. In 
b) the universal TA curve at high Ti , T2 was subtracted to 
remove completely the TA background. This allows clear ob- 
servation of evolution of the sum-gap peak and the double-gap 
dip well above the phase-coherent t^^°-'"^. 



mesas are represented by symbols of the same type and 
color. 

Dashed-dotted lines represent the TA barrier, ob- 
tained from the zero-bias resistance using Eq.([2]) Uta = 
kBTln{Ro/T). It is seen that Uta is practically constant 
at T > Tc. The sudden fall of Uta marks the supercon- 
ducting transition. 

Crosses in Figs. 10-12a) represent crossing points 'io], 
marked in Figs. 2 and 4b). In agreement with Fig. 6 
they occur at eV/N ^ 2Uta- 

Open symbols in Figs. 9-11 a) show the hump voltage 
[2^ . It also represents the TA barrier and is roughly 



T— independent at T > (26l. \2T\. The hump appears 
at slightly higher voltage than 2Uta/g, also in agreement 
with Fig. 6. 

From Figs. 10-12 a) it is seen that in most underdoped 
crystals a slight deviation of Uta downwards, and the 
crossing point upwards occur below some temperature 
T* ^ 130— IbQK. Dashed and dotted vertical lines mark 
the Tf'"'" and T* , respectively. 

Solid symbols in Figs. 9-12 a) represent the main 
experimental result of this work: T— dependencies of 
sum-gap voltages eVgg/N = 2 A. Data points were ob- 
tained from the peak and half the double-gap dip in 
dl/dV (larger symbols), or T— differential characteris- 
tics (smaller symbols). The latter allows us to trace the 
gap up to considerably higher temperatures, than before. 

In agreement with previous reports [2^ [27j the bulk 
gap considerably decreases at T ^ Tc for all doping 
levels. Simultaneously, we observe that the hump also 
becomes T— dependent at T < Tc [13]. However, it 
moves approximately two times slower than the sum-gap 
peak and is approximately described by the expression 
eVhunipiT < Tc) ~ 2Uta + A(T). This indicates that 
the c— axis TA barrier remains intact by the supercon- 
ducting transition [26] and continue to hinder the QP 
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FIG. 9: (Color online), a) Temperature dependence of char- 
acteristic spectroscopic features for three mesas on the near 
optimally doped SMa crystal. Data for the same mesas are 
shown by symbols of the same type and color. Open symbols 
- the TA hump, the dashed-dotted curve - 2Uta- Downturn 
of Uta marks the c— axis resistive transition, TP'^'^'"' at which 
c— axis phase coherence is established is marked by the dashed 
vertical line. Solid symbols - the sum-gap peak, it is seen that 
the measured gap is distorted with increasing mesa area ex- 
actly as expected for self-heating. Inset in a) shows power 
at the peak for the same mesas, b) The heat-compensated, 
T— dependence of the bulk energy gap. Solid line shows that 
it is very well described by the conventional BCS temper- 
ature dependence with the mean-field critical temperature 
T^^ = 96K. No pseudogap is observed above T^^ . Inset 
in b) shows the effective thermal resistance for the two larger 
mesas. 



transport at T < Tc 



A. Size dependence 

Size-dependence of ITS unambiguously reveals the ex- 
tent of self-heat distortion [5^ . Such data is presented in 
Figs. 9 and 12. Dissipation powers at the sum-gap peak 
Ppeak for all studied mesas are shown in insets of Figs. 9 
a) and 10-12 b). Ppeak scales with the mesa area. 

Solid symbols in Fig. 9 a) show measured Vsg for three 



mesas on the SMa crystal. Size dependence of ITS for 
this crystal was reported in Ref.fS^]. It was shown that 
for mesas with A < IS/im^ the measured gap becomes 
size- independent. Therefore, data for the smallest mesa 
with A ~ 3.6/im^ represents the genuine, undistorted 
bulk gap A(r), as concluded in sec. Ill A. The solid 
line in Fig. 9 a) shows that it is very well described 
by the conventional mean-field BCS temperature depen- 
dence Abcs{T). The same is true for the small mesa on 
the S92 crystal in Fig. 10a), which was also identified as 
heating-free in sec. Ill A. 

Let's now consider the two larger mesas in Fig. 9a). It 
is seen that measured gaps become progressively smaller 
with increasing A and Ppeak- The observed deviation 
from the genuine A(T) is perfectly consistent with self- 
heat distortion, as shown in Fig. 15 d). The temperature 
rise 5T is equal to the horizontal shift of the measured 
gap with respect to the undistorted A(T), as indicated 
in Fig. 9 a). 

Now we can directly calculate the effective thermal re- 
sistance of the mesas: RrhiT) = ST/Ppeak{T). The ob- 
tained RTh{T) appeared to be approximately the same 
for both mesas and is shown in the inset of Fig. 9 b). The 
values of Rth are ranging from ^ 25K/mW at 4:.2K to 
^ lOK/mW at Tc, consistent with direct in-situ mea- 
surements in Ref. [sj. Note that in Ref. [5^ a separate 
thermometer was employed for measuring the mesa tem- 
perature. Therefore, agreement in obtained RTh{T) in 
both cases indicates that there is no major thermal gra- 
dient along the crystal near the mesa. This is consistent 
with the conclusion of Ref. [HOl that heat transport from 
the mesa is dominated by ballistic flow of non-equilibrium 
phonons. In this case the actual heating starts only deep 
inside the crystal, phononic mean-free path below the 
mesa, making both heating and in-plane thermal gradi- 
ent in the mesa small. 



B. Self- heat compensation 

The knowledge of the thermal resistance allow us to 
recover the genuine temperature dependence of the gap 
even for moderately large mesas. In Fig. 9 b) gap 
values for all three mesas from SMa are plotted as a 
function of the actual mesa temperature Teorr ^ T + 
Ppeak {T)RTh{T), with Ppeak{T) and RxhiT) from insets 
in Figs. 9 a) and b), respectively. It is seen that the gen- 
uine T— dependence of the gap is recovered also for larger 
mesas after such self-heat compensation. The solid line 
in Fig. 9 b) shows that the genuine A(T) is perfectly de- 
scribed by the BCS dependence with T^^ = 96K, which 
i approximately equal to the optimal Tc of our Bi(Y)-2212 
crystals. It is close to J'°"'"=* ~ 95K and slightly higher 
than the phase-coherent Tc ~ 91.5K for this crystal, as 
shown in Fig. 13. 

We can also check the self-consistency of the conclusion 
that ITS characteristics of the smallest mesas on SMa 
and S92 are not distorted by self-heating. As seen from 
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FIG. 10: (Color online). The same as in Fig. 9 for the small 
mesa on the moderately under S92 crystal. Crosses in a) 
represent the crossing point of dl /dViV) curves at two con- 
secutive temperatures, see Fig. 4 b). Filled triangles rep- 
resent half the double-gap dip. It is seen that signatures of 
the energy gap (the pseudogap) survive up to T* ~ HQK (see 
Fig. 7). Solid and dashed lines in b) show that A(T) is equally 
well described by the pure BCS-dependence or the combined 
gap between the BCS gap and T— independent pseudogap, 
shown by the dashed horizontal line. However, in both case 
the mean-field superconducting temperature is smaller than 
that for the optimally doped crystal in Fig. 9. 

insets in Figs. 9 a) and 10b), the maximum Ppeak for 
those mesas is ~ 0.2mW at the lowest T. Therefore, 
at the lowest T the maximum self-heating 6T does not 
exceed few K, which does not affect the measured gap. 
Ppeak rapidly decreases with increasing T. At Tc/2, Ppeak 
reduces by half to ^ O.lmW^ and 6T ^ 2K. Close to 
Tc, self-heating becomes negligible even for moderately 
large mesas. For the smallest mesas from Figs. 9 and 10, 
PpeakiT ~ Tc) ~ m^iW and 5T is sub-Kelvin. Therefore, 
the measured gaps for the two smallest mesas are indeed 
undistorted by self-heating. 

Figs. 10-12 b) represent self-heat compensated A(T) 
for moderately underdoped mesas. The same RxhiT) 
was used for self-heat compensation. In all cases, the 
recovered A(T) can be very well described by the mean- 
field BCS dependence (solid lines) with T"'-^ that is larger 



FIG. 11: (Color online). The same as in Fig. 10 for the S43 
crystal. It is seen that the superconducting transition at Tc 
even in moderately underdoped HTSC has a mean-field char- 
acter and is accompanied by opening of the superconducting 
gap. 

than Tf'^"-^'^ and close to T""**^*. Those values are sum- 
marized in Table-I. 



C. Phase-coherence in c— axis and a6— plane 

In recent years there were many speculations about 
persistence of phase-incoherent superconductivity up to 
very high temperatures above Tc in underdoped HTSC. 
Therefore, it is necessary to clarify the definition of Tc 
and the difference between c— axis and a6— plane phase 
coherence. 

The c— axis phase coherence is caused by weak Joseph- 
son coupling between CuO planes: Ej = (ft/2e)/c, where 
Ic is the total Josephson critical current of the junction. 
The Ic and the Josephson coupling can be easily sup- 
pressed by small magnetic fields, which do not affect 
superconductivity of planes. Moreover, the Josephson 
coupling can be simply reduced by reducing the junction 
area. At T > Tj = Ej/ks 23.8K{Ic/ fiA), thermal 
fluctuations destroy phase coherence and the junction en- 
ter in the phase diffusion state with non-zero resistance at 
zero current [13, HH . Suppression of the critical current 
by thermal fluctuations at T ^ Tc was indeed observed 
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FIG. 12: (Color online). The same as in Fig.lO for the S82 
crystal. Data for mesas with different area at T ~ 6K are 
shown. Gray area in b) highlights the T— region without 
c— axis phase coherence but with persisting gap and in-plane 
superconductivity. 



in small mesas [13|, |6l|. Therefore, the c— axis phase co- 
herence in small Bi-2212 mesas is not a material property 
and should not be confused with a much more robust in- 
plane phase coherence. 

To clarify the difference between in-plane and c— axis 
phase coherence, in Fig. 13 a) I show details of in- 
plane and c— axis resistive transitions, measured on a 
slightly underdoped crystal SIz, similar to SMa. It is 
seen that the c— axis phase coherence is established at 
rpphase ^ gg,^K. The onset of c— axis transition occurs at 
T°'^^'^*(c—axis) ~ 93iir, which is close to the middle point 
of the in-plane transition. The onset of a6— plane resistive 
transition occurs at T°"''^*{ab) ~ 96K, which coincides 
with the mean-field critical temperature TJ"-^ ~ 96K, ob- 
tained by extrapolation of A(r) using BCS temperature 
dependence, as shown in Fig. 13 b) 



D. The pseudogap 

The main difference between near optimally doped 
(Fig. 9) and moderately underdoped (Figs. 10-12) crys- 
tals is the persistence of some signature of the resid- 



FIG. 13: (Color online), a) In-plane Rab and c— a^cis Rc re- 
sistive transitions in a slightly underdoped crystal, similar to 
SMa. Panel b) reproduces A(r) from Fig. 9 b) along with the 
resistive transition for a mesa on SMa. It is seen that the onset 
of in-plane superconductivity coinsides with the extrapolated 
mean-field critical temperature, the onset of c— axis transi- 
tion is close to the middle of the in-plane resistive transition, 
while the temperature at which c— axis phase coherence is es- 
tablished is further reduced due to thermal fluctuations. Gray 
area in b) marks the superconducting region without phase 
coherence in the c— axis direction. 



ual energy gap (the pseudogap Apc) at T > Tc up to 
T* - 130 - 150K (see Figs. 7, 8). In the same range 
Tc < T < T* we observe an upturn of the crossing voltage 
and downturn of Uta obtained from zero bias resistance. 
Thus the pseudogap is seen at all bias levels. 

From Figs. 10-12 b) it is seen that the PG is almost 
T— independent near T™/ in contrast to the strong A(r) 
dependence on the superconducting side. It is also seen 
that the PG is merging with the superconducting gap at 
Tc- The overall T— dependence of the gap both above 
and below Tc is well described by the combined gap ex- 
pression: 



Aco,nb{T) = ^^BCS{TY + A|g, (4) 
with constant Apg. Fits of Eq.(I3|) to experimental 
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A(T) are shown by dashed blue hnes in Figs. 10-12 b). 
Below Tc they are equally good as BCS fits without the 
PG, but assume lower T^-^, indicated in the figures and 
Table-I. 

Can this pseudogap be related to superconductivity? 
The answer is rather straightforward. The very fact that 
it forms the combined gap with the superconducting gap 
clearly indicates that it represents another order param- 
eter, co-existing and competing with superconductivity. 

Our present data, together with earlier observation of 
coexistence of the superconducting gap and the pseudo- 
gap ai T < Tn [2^, are consistent with recent ARPES 
experiments [2], which demonstrated that the supercon- 
ducting gap along the Fermi arc closes at at all doping 
levels, while the gap along antinodal directions (the pseu- 
dogap) remains relatively T— independent at T^. 



V. DISCUSSION AND CONCLUSIONS 

I have shown that technical problems of Intrinsic Tun- 
neling Spectroscopy, such as self-heating (see Appendix- 
A) and in-plane resistance (see Appendix-B), can be ef- 
fectively obviated by reducing mesa size. It is instructive 
to remind that exactly the same technical problems were 
encountered at the early stage of experimental studies of 
conventional LTSC tunnel junctions fssj and were also 
solved by junctions miniaturization, which reduces the 
area-to-perimeter ratio of the junctions. 

I have also shown that self-heating is a trivial phe- 
nomenon. It can be easily accounted for and, in case 
of moderate heating, compensated, so that the genuine 
A(r) can be recovered. A rule of thumb for "moderate" 
self-heating is that it should not be obvious in the IVCs: 
the QP branches should reman periodic and there should 
be no back-bending at Vgg ■ All mesas studied here fall in 
this category. 

The main result of this work is the uncovered gen- 
uine temperature dependence of the bulk energy gap in 
Bi-2212. For all studied doping levels, A(r) exhibits a 
strong T— dependence, and the superconducting part of 
it unmistakably tends to vanish in the mean-field BCS 
manner. For slightly overdoped Bi(Y)-2212 crystals from 
the same batch this was shown in Ref.|,27j. Here I have 
focused on the underdoped side and have shown that the 
gap vanishes exactly in the BCS manner also in slightly 
underdoped (Fig. 9b) and moderately underdoped (Figs. 
10-12) mesas. 

Those results are strikingly different from the complete 
T— independent surface gap, reported in STM experi- 
ments [5]. I want to emphasize that this discrepancy 
can not be attributed to self-heating and must find an- 
other explanation. Indeed, numerical simulations in Fig. 
15 clearly show that the trivial self-heating simply can 
not "hide" the qualitative A(r) dependence. For exam- 
ple, there is no way in which one can get the vanishing 
"measured" (self-heating affected) gap if the true gap is 
T— independent. Furthermore, self-heating becomes in- 



significant at elevated T because Ppeak{T — > T^f) — > 0. 
On the other hand, STM characteristics behave similar 
to the T— independent hump feature in ITS, which is the 
consequence of the c— axis thermal activation barrier. As 
discussed in sec. III-B, Uta is most likely the property of 
the blocking BiO plane, which is probed by STM, rather 
than superconducting CuO planes. The dramatic dif- 
ference between STM spectra on BiO and CuO surfaces 
was indeed reported [3]. This highlights the significance 
of spectroscopic information from bulk CuO planes, ob- 
tained here. 

Main conclusions, that can be drawn from ITS data 
are: 

(i) The superconducting transition, even at moderate 
underdoping, is due to conventional mean-field phase 
transition rather than destruction of phase coherence 
without amplitude fluctuations. This confronts specula- 
tions about persistence of the "precursor" superconduct- 
ing state in the extended T— region above Tc- 

(ii) The mean-field superconducting critical tempera- 
ture decreases with underdoping. Thus, HTSC does not 
become stronger with approaching the undoped antifer- 
romagnetic state. 

(iii) The pseudogap co-exists and competes with su- 
perconductivity and disappears at optimal doping. 

An important consequence of those conclusions is that 
high temperature superconductivity is strongest at op- 
timal doping and becoming weaker with underdoping. 
This is consistent with the decrease of the upper criti- 
cal field [2§| with underdoping. Therefore, our observa- 
tions support the idea that the mechanism of HTSC is 
intimately connected to a Quantum Critical Point near 
optimal doping [3, Ig^ . rather than closeness to the an- 
tiferromagnetic state |30| . 
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APPENDIX A: ANALYSIS OF SELF-HEATING 
IN ITS 

Despite relative simplicity of self-heating phenomenon 
(it is certainly the most trivial problem in HTSC spec- 
troscopy), discussion of self-heating in ITS has caused a 
considerable confusion, a large part of which has been 
caused by a series of pubhcations by V.Zavaritsky [63j . 
in which he "explained" the non-linearity of ITS charac- 
teristics by assuming that there is no intrinsic Josephson 
effect. Irrelevance of that model for ITS was discussed in 
Refs. [511. 



A certain confusion might be also caused by a large 
spread in thermal resistances, Rrh: reported by differ- 
ent groups [13, [El, [il, EH]. For the sake of clarity 
it should be emphasized that those measurements were 
made on samples of different geometries. It is clear that 
Rxh depends strongly on the geometry [H, H^], and is 
much larger in suspended junctions with poor thermal 
link to the substrate [6^ than in the case when both 
top and bottom surfaces of the junctions are well ther- 
mally anchored to the heat bath [13] ■ For mesa struc- 
tures similar to those used in this study (a few /im in- 
plane size, containing ~ 10 IJJ), there is a consensus 
that RTh{4:.2K) ~ 30 - 70K/mW (depending on bias) 
and RTh{90K) ~ 5 - lOK/mW 0. Larger val- 
ues Rxh > lOOK / mW claimed by some authors [5a| are 
unrealistic for our mesas because they can withstand dis- 
sipated powers in excess of lOmW without being melted. 

Yet, talking about a typical value of Rxh is equally 
senseless as talking about a typical value of a contact 
(Maxwell) electrical resistance: both depend on the ge- 
ometry. Therefore, reduction of mesa siz es p rovides a 
simple way for reduction of self- heating [5^. Conse- 
quently, variation of dl /dV characteristics with the junc- 
tion size and geometry provides an unambiguous way of 
discriminating artifacts of self-heating from the spectro- 
scopic features [50]. 



1. Peak splitting in non- uniform junctions 

As discussed in Ref.]51], atomic separation between 
IJ Js precludes any substantial temperature difference be- 
tween them. Thus, all junctions in a mesa warm up syn- 
chronously and there may be only one collective artefact 
of heating for all IJJs in the mesa. To the contrary, if 
the peak is the sum-gap singularity - it is an individual 
property of each IJJ. If junctions are not perfectly iden- 
tical, the peak in dl/dV will split in up to N sub-peaks. 
Indeed, peak splitting is quite often observed in exper- 
iment and was already reported in Ref.[66] along with 
supporting numerical simulations. 

In Fig. 14 the IVC at T = 30K from Fig.2 is re-plotted 
in a semi-logarithmic scale V vs In/. As discussed in 
Ref.fi^ the IVC becomes nearly linear on such a scale. 
Thin lines are multiple integers of the last branch di- 
vided by iV = 9. Coincidence of those with quasiparticle 
branches indicates good periodicity of the latter. How- 
ever, a minor non-uniformity is observed as a gradual 
increase of Ic with the brunch number. The total spread 
of the critical current Ale from the second to the last IJJ 
is marked in the Figure. The red line shows dI/dV{I) 
for the same IVC. A small splitting of the peak A/^g is 
seen. Thus, the peak is not a collective phenomenon of 
the whole mesa, but is a genuine characteristic of each 
individual IJJs. From Fig. 14 it is seen that Algg has 
approximately the same width in the logarithmic scale 
as Ale- Therefore the splitting is proportional to the 
difference in critical currents of IJJs in the mesa and is 




/(mA) 

FIG. 14: (Color online). The IVC at T = 30K from Fig.2 in 
a semi-log scale. Thin lines are multiple integers, indicating 
good periodicity of quasiparticle branches. However, a small 
nonuniformity of the junctions is seen from the spread of criti- 
cal currents Ale, which leads to splitting A Jag of the sum-gap 
peak in dl /dV{I) (red line) because different junctions reach 
Vag at slightly different currents. 

due to the corresponding spread in currents I^g at which 
individual junctions reach Vsg. 

Thus, non-uniformity of junctions, although usually 
unwanted for ITS, helps to understand the origin of the 
peak in ITS characteristics of small Bi-2212 mesas. 



2. How self-heating affects I — V characteristics. 

How self- heating can distort the IVC's of Josephson 
junctions is obvious: since self-heating rises the effective 
T it may affect the IVC only via T— dependent parame- 
ters. There are three such parameters: 

i) the quasiparticle resistance, 

ii) the superconducting switching current, 

iii) the superconducting gap. 

They will affect the IVC of mesas, containing several 
stacked IJJs, in the following manner: 

The consecutive increase of T upon sequential switch- 
ing of IJJs from the superconducting to the resistive state 
will distort the periodicity of quasiparticle branches. 
Each consecutive QP branch will have a smaller QP re- 
sistance (smaller V at given /) and smaller switching 
current. This type of distortion becomes clearly visible 
(at base T = 4:.2K) when AT > 20K 

For better understanding of the influence of self- 
heating on IVCs of Josephson junctions, in Fig. 15 I 
reproduce the results of numerical simulation of such the 
distortion, made specifically for the case of Bi-2212 mesa 
with the corresponding T— dependent parameters (see 
Ref.[6^ for details). Fig. 15 a) shows a set of undistorted 
IVC's at different T for coherent, directional, d— wave 
tunneling with some trial A(T), shown by the solid line 
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FIG. 15: (Color online). Simulated distortion of SIS tun- 
neling characteristics by strong self-heating, a) Undistorted 
IVC's at different T . Simulations were made for typical pa- 
rameters of our mesas, T— dependent thermal conductivity of 
Bi-2212, and for coherent, directional, d— wave tunneling, b) 
distorted IVC's at the same base T; c) The mesa tempera- 
ture as a function of bias, d) T— dependence of the genuine 
superconducting gap (solid line) and the "measured" gap ob- 
tained from distorted IVC's (dashed curve). Note that even 
strong self-heating (T reaches rc/2 at Kg at 4.2K) does not 
cause considerable distortion of the measured gap. Data from 
Ref.lii 



in Fig. 15 d). Panels b) and c) show the distorted IVCs 
and the actual junction temperature, respectively. It is 
seen that combination of self-heating and T— dependence 
of A may lead to appearance of back-bending of the IVC 
at the sum-gap knee. The dashed line in panel d) repre- 
sents the "measured" gap obtained from the peak in dis- 
torted dl/dV characteristics. Remarkably, the deviation 
from the true A(T) is marginal, despite large self-heating, 
AT ~ Tc/2 at A.2K\ Numerical simulations has shown 
that even self-heating up to Tc at the sum-gap knee does 
not cause principle changes in the behavior of the "mea- 
sured" gap. The robustness of the measured gap with 
respect to self-heating is due to the flat T— dependence 
of the superconducting gap at T < Tc/2 and to simulta- 
neous vanishing of dissipation power at Vgg together with 
A(T) at T ^ Tc, as shown in insets of Figs. 9-12. 



3. Heating or non-equilibrium phenomena? 

Finally, it is important to emphasize that the concept 
of heat diffusion is inapplicable for small Bi-2212 mesas 
containing only few atomic layers. The phonon transport 
in this case is ballistic [13, [63| and the energy flow from 
the mesa is determined not by collisions between the tun- 
neled non-equilibrium QP with thermal phonons, but by 



spontaneous emission of a phonon upon relaxation of the 
non-equilibrium QP ^38j . This process is not hindered at 
T = 0. Therefore, the effective Rrh (and self-heating) 
can be much smaller because it is not limited by poor 
thermal conductivity at T = 0, but is determined by the 
fast, almost T— independent, non-equilibrium QP relax- 
ation time. The concept of self-heating becomes adequate 
only in the bulk of the Bi-2212 crystal, where the dissi- 
pation power density and the temperature rise are much 
smaller due to the much larger area of the crystal. For 
more details see the discussion in Ref. [s^l. The non- 
equilibrium energy transfer channel is specific for atomic 
scale intrinsic Josephson junctions made of perfect single 
crystals. It can explain a remarkably low self-heating at 
very high bias [s^ . 



APPENDIX B: ARTIFACTS OF IN-PLANE 
RESISTANCE 

In experiments on large mesas [s^ or suspended struc- 
tures [65| no clear Ohmic tunneling resistance could be 
observed. Instead a negative differential resistance (acute 
back-bending) progresses at high bias. Similar behavior 
is also observed on moderately large mesas, when mea- 
surements are made in the four-probe configuration [68| . 
This is in stark contrast to three-probe measurements 
on small mesas, reported here, see Fig. 2 and Refs. 
[2i,|2l|4i,[53. 

The continuous negative resistance is not described by 
self-heating in tunnel junctions because those should al- 
ways reach the T— independent positive tunnel resistance 
at high bias. Thus self-heating is not the primary cause 
of the acute back-bending. Rather, the negative differ- 
ential resistance, observed in large structures is caused 
by the loss of equipotentiality of CuO planes upon which 
the measurement of IVCs in mesa structures is relying. 
The latter can be triggered by development of a hot spot 
(stI . but can even occur without self- heating due to in- 
plane resistive transition of CuO planes when the applied 
current exceeds the critical current of the CuO plane [11] . 



1. Acute back-bending without heating and the 
difference between 3-and 4-probe measurements 

Here I demonstrate how acute, and not recovering 
back-bending develops in the IVC of the mesa as a re- 
sult of the finite in-plane resistance and without any 
self-heating. I also explain the difference between four- 
probe measurements, that exhibits acute back-bending 
[6a |. and three-probe measurements - that don't. 

The top panel in Fig. 16 show the sketch of the mesa 
structure. The current / is biased through the top elec- 
trode and returned through the crystal, which serves as 
the bottom electrode. The voltage can be measured ei- 
ther in four-probe V/^p or thee-probe Vijp configuration, 
as shown in the Figure. 
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2 

V /[/*/?] 



v.. = 



IR1R2 

{Ri + i?2)a 



,t + t+2cosha 
sinh a 



(Bl) 



where a = 4/ ^^"^^^ . In case i?i = i?2 it coincides with 



the result of Ref.[^. 

In the three-probe configuration the measured voltage 
is Vsp = Vi{0) - ¥2(0), which yields: 



i?2 + Ri cosh a 
sinh a 



(B2) 



The main panel in Fig. 16 shows calculated four-probe 
(solid lines) and thrcc-probe (dashed lines) IVCs for the 
case of identical superconducting electrodes Ri = R2 = 
Rab{I) with the in-plane critical current and the IVC 
as shown in the inset. For Rab = the measured IVC 
coincides with real c— axis IVC of the mesa. But sub- 
stantial deviations occur when the bias current exceeds 
1°:^ and electrodes become resistive. 



Limitations on the mesa size 



FIG. 16: (Color online). Top panel shows a sketch of current 
distribution in case of resistive electrodes and contact config- 
uration for three- and four-probe measurements. The main 
panel shows simulated I — V curves in the three-probe (dashed 
lines) and four-probe (solid lines) configurations for different 
in-plane electrode resistances. It is seen that in-plane resis- 
tance distorts the measurement of junction characteristics nd 
that acute back-bending can develope in the four-probe case 
for high in-plane resistance. Inset show the in-plane IVC. 

If electrodes (both top and bottom) are resistive, then 
there is a voltage gradient along the electrodes and 
the bias current is distributed non-uniformly within the 
mesa, as shown by thin vertical arrows. Let Ri^2 be the 
in-plane resistance of top and bottom electrodes within 
the mesa area > a: > a, and R^ the c— axis resistance 
of the mesa itself. Following Ref.jsl], the current and 
voltage distribution in this case can be described by the 
system of equations: 

dx a ' 

dV2 R2 , 
—r- 

dx a 

h{x)+l2{x)^I, 

Vi{x)-V2{x) = -Rca^, 
dx 

with the boundary conditions: /i(0) = / and Ii{a) = 

0. 

In the four-probe configuration, the measured voltage 
is V4p = Vi{a) — ¥2(0) , which yields: 



I want to emphasize that large mesas are not suitable 
for ITS even in the absence of self-heating. To probe the 
gap, one should be able to reach the sum-gap voltage, 
Vsg ~ 60meV per IJJ, without loosing the equipotential- 
ity of the bottom CuO plane, which is used as the return 
current lead and the second voltage electrode. 

Let's estimate the maximum mesa size, suitable for 
ITS, in the absence of self- heating. Consider a square 
mesa with the in-plane size a. The bias current, required 
for reaching the sum-gap voltage is Igg — Vsg/{pcs/oP'), 
where pc — SOSlcm is the c— axis tunnel (large bias) re- 
sistivity, and s ~ 1.5nm is the interlayer spacing. This 
current is flowing through the perimeter of the last IJJ 
into the bottom CuO layer and should not exceed the in- 
plane critical current of that layer. Provided the in-plane 
critical current density is J"** ~ lO^A/cvn? [6^, the in 
plane critical current of the bottom CuO plane through 
the perimeter of the mesa is ^ Ufsa. Therefore, the 
mesa size should be smaller than amax — ^Jc'^ Pc/Vsg — 
4.5/im. Thus, miniaturization is essential for ITS. Oth- 
erwise the return current and voltage contacts are no 
longer equipotential, leading to the negative measured 
differential resistance. 

To cause a substantial distortion, the total in-plane re- 
sistance should be larger than the mesa resistance. Let's 
see if this is the case for Bi-2212. The in-plane resistance 
of the bottom CuO plane in the square a x a is: 

Rab{^) = Paba/sa - 10-'^(ncm)/1.5(nm) = 667(f7). 

The c— axis mesa resistance is Rc = PcNs/a? ~ 
30(r2cm)1.5(nTO)7V/a2 = AhQN / [a{pm)f {n) . 

For a mesa 5 x bpw? with iV = 10 IJJs, Rc = ISOfi is 
about four times smaller than i?afc(n). Thus we see that 
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distortion by in-plane resistance can indeed be significant 
for larger mesas with a small number of junctions. 

The proposed model explains why four- and three- 
probe measurements of Bi-2212 mesas may be very dif- 
ferent. From Fig. 16 it is seen that four- and three-probe 
measurements react differently on in-plane resistivity. 
For the three-probe case, it just leads to appearance of 
an additional series resistance. But for the four-probe 
configuration it may lead to development of the acute 



back-bending, precluding any spectroscopic analysis. 

The model can also explain a strangebehavior of mesas 
with very small amount of junctions [69|. From the esti- 
mations above it follows that for a mesa 5 x with 
only one IJJ the in-plane resistance is about 40 times 
larger than the mesa resistance, which make such mesas 
extremely prone to distortion by in-plane resistance and, 
probably, not suitable for ITS. 
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